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. ABSTRACT 

<D 

pL^ ' Context. The formation mechanism of hot horizontal branch (HB) stars is still one of the most uncertain points of stellar evolution 

theories. In the past decade, models based on their binary origin have been very successful in reproducing the properties of field 
subdwarf-B stars, but the observations of their analogues in globular clusters has posed new problems. In addition, the discovery of 
multiple populations offered an appealing alternative scenario for the formation of these stars. 

Aims. We search for binaries of period p <200 days among a sample of 83 blue horizontal branch stars (T c ff= 12 000-22 000 K) in 
NGC 2808, a cluster known to host three distinct stellar populations and a multimodal horizontal branch. The final sample, after the 
rejection of stars with incomplete observations or poor quality data, consists of 64 targets. 
Methods. The radial velocity of the targets was measured in fourteen epochs, spanning a temporal interval of ~75 days. The significant 
variations were identified by means of a detailed error analysis and a statistical study. 

Results. We detect no RV variable object among stars cooler than the photometric Gl gap at ~17000 K, while two close (p <10 
days) and two intermediate-period (p= 10-50 days) systems are found among hotter targets. The close and intermediate-period binary 
fraction for stars cooler than the gap are / c <5% and f ip <10%, respectively, with 95% confidence. The most probable values among 
hotter stars are / c ~20% and /| p ~30%, but the 90%- confidence level intervals are still large (6-42% and 1 1-72%, respectively). 
Conclusions. The Gl gap appears as a discontinuity in the binary faction along the HB, with a higher incidence of binaries among 
hotter stars, but a constant increase in / with temperature rather than a discontinuity cannot be excluded from our observations. We 
also find that intermediate-period binaries, never investigated before among cluster HB stars, could play an important role among 
hotter stars, being more than ~ 15-20% of the hottest stars of our sample. Our results, compared with previous estimates for other 
clusters, indicate that f c among hot HB stars is most probably higher for younger clusters, confirming the recently proposed age- 
f c relation. However, the large observed difference in binary fraction between clusters (e.g. NGC 2808 and NGC 6752) is still not 
reproduced by binary population synthesis models. 
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1. Introduction NGC 6752. lMomanvetal]|2q02l). while in oth ers it is clearly 

TT . „ , /TT „, • ^, • , , multimodal (e.g. NGC2808. ISosin et al"1ll997l) . This observa- 

Honzontal branch (HB) stars in Galactic globular clusters (GCs) tional icture &m lacks m comprehensiori5 as a conseqU ence 

are old post- He flash stars of low initial mass (0.7-0.9 M ) that, of Qur unders tanding of the formation mechanism of HB 

after the exhausnon of hydrogen in the steflar core and their as- starg fa GCfj The HB morphol has been linkedi a 

jH . cension along the red giant ^branch, eventually ignited helium tQ duster age (DoUer gt aJ 2010)> duster concentration 

C?_. d Hoyle & Schwarzschild| | l955HFaulkneHl966| ). (Rlsi pecd £t aJ i 993) , ste llar rotation (PetersolE^, duster 

GCs display large differences in the HB morphology (see, for mass ( Recio . Blanco et al . 20 06), helium, and the environment of 

example,| Piottoetal|2002|) . The first parameter responsible for formation (Fraix . B urnet et al. 2009), but none of the proposed 

this phenomenon is metalhcity, but it alone cannot account for secQnd parameters cou i d satisfa ctorily reproduce the complex 

the complex observational picture (the so-called "second param- observed beh avior (see Catelan 200l for a review). In particu- 

eter problem , [Sandage & Wildey || 1967t | van den Bergh| | 1967| ). ^ ±e mQSt challenging task is t0 account for the formation of 

While some clusters contain only red HB stars cooler than extreme horizontal branch (EHB ) stars at the faint hotter end of 

the RR-Lyrae gap, others host a large population of blue He- HBfj (Teff > 20000 K)> observed even in hi h meta iii city dus . 

burning stars extending even beyond the canonical end of the tgrsfgjjR ich et Zl[l997h and old open clusters (NGC 6791, 

HB at -35 000 K (e.g. | Moehler et alj |2004D. Even more puz- g USO n et al. 2006). Stars hotter than this critical temperature do 

zling, in the color-magnitude diagram (CMD) of some clus- nQt haye an extemal envel massive enou h to sustain the 

ters the HB appears continuous in its whole extension (e.g. ^ H . burning; and after the exhaustion of helium in the core 

I _ T ; ; ., , r „„ ,, . they evolve directly to the white dwarf (WD) cooling sequence, 

Based on observations with the ESO Very Large Telescope at . , ,. , ^t-, ' 

Paranal Observatory, Chile (proposal ID 078.D-0825) wlthout ascending the asymptotic giant branch (AGB manque 
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stars. iGre ggio & Renzini 1990). They are extensively observed 
and studied in the Galacti c field, identified as the so-called sub- 
dwarf B-type (sdB) stars (lGreensteinlll971t ICaloilll972l iHeberl 
119861) . However, in GCs they are still poorly studied because of 
their faintness , and many questions still a wait an answer (see 
ICatelari 120091: iMoni Bidin & Piottd [20ToL for recent reviews). 
Many single-star evolutionary channels have been invoked to ex- 
plain EHB sta r formation in GCs, i ncluding interactio ns with a 
close planet dSokerfl!998l see also ISilvotti etaTll2007l). helium 
mixing driven by eith er internal rotation (ISweigart & Men gel 
119791; IS weigarl 1 1997b or stellar encounters dSuda et all 120071) . 
and close encou nters with a central intermediate-mass black hole 
(lMiocchill2007h . 

The dynamical interactions inside binary systems were pro- 
posed early on to be r esponsible for the heavy mass-loss required 
to for m a EHB star (Menge fet al.lll976t iTutukov & Iungelsonl 
1 19871) . In the past decade, the "binary scenario" has achieved 
many observational and theoretical successes among field sdB 
stars, and is now widely accepted as the most satisfactory ex- 
planation of their formation. The bin ary population synthesis 
model o flHan etali (2002. 20031 120071) could reproduce their ob- 
servational properties in great detail, although a small fraction 
of prog enies of single st ars is probably required for a perfect 
match (Lisker et al. 2005). Han's model considers three main 



formation channels: the stable Roche Lobe Overflow (RLOF), 
which produces sdB's in wide binaries; the common envelope 
(CE) channel, which forms close systems; and the merging 
of two WDs, whose progenies are single stars. On the other 
hand, many surveys have con firmed that a large fraction of field 
sdB stars reside in binaries (Ferguson et al.lll984l lAllard et al.l 
1991; lUlla & Theilll [1991 lAznar Cuadrado & Jeffervl 120011; 
Maxted et alj20qUl Williams et alj200UlReed & Stiening|2 004: 
Napiwotzki et al.l 120041) . and sdB's in close syste ms with pe- 



riods shorter than ten days are very common (|Moran et aT 
1999tlSaffer et al.lll998tlHeber et aUl2002l:lMorales-Rueda et al 



2003), although the exact cl ose-binary fraction is st ill uncer- 
tain, ranging from 40-45% dNapiwotzki et al.l |2004 to 70% 
dMaxtedetalJl2001l) . 

Observations of EHB stars in GCs have so far presented 
a challenge to the binary scenario, at variance with its well- 
established successes for field stars. The first surveys sur- 
prisingly revealed a lack of EHB close systems in GCs 
dMoni Bidin et al.ll2006ir2009bl) . lMoni Bidin et alj J2008a) fixed 
at 4% the most probable value of the EHB close-binary fraction 
(/ c ) in NGC 6752, proposing that a decrease in / c with the age 
of the stellar population should be a natural expec t ation of the 
binary scenario. The detailed calculations of Hani (j2008) con- 
firmed that the CE channel becomes very inefficient after the 
first few Gyrs, and the WD-WD merging should be the predom- 
inant mechanism for the formation of EHB stars in old stellar 
systems. As a consequence, cluster EHBs should be principally 
single-star pro ducts of merging, and clo se systems should be 
rare. However, Moni Bid in et al . (2009a) measured a higher f c 
in M80 and NGC 5986 (12% a nd 25%, res pectively), and ar- 
gued that the models proposed bv lHanl (120081) cannot simultane- 
ously account for the low f c measured in NGC 6752 and these 
much higher values found in clusters only 1-2 Gyr younger. 
Unfortunately, their results suffered too large uncertainties to be 
conclusive. In summary, the binary scenario has not been dis- 
proved, but its ability to reproduce all the observations in GCs 
still has to be demonstrated. 

An alternative model of EHB star formation has received 
great attention in recent years: one incorporating the primordial 
helium enhancement. A super-solar surface helium abundance 
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Fig. 1. Position of the targets in the cluster color-magnitude di- 
agrams. Upper panel: V vs (U - V) diagram, where the arrows 
indicate the location of the two gaps discussed in the text (Gl 
and G2), and the magnitude corresponding to T=20 000 K (the 
formal definition of the EHB). Lower panel: V vs(B-V) d i agram 
with superimposed the ZAHB model from ICassisi et al.l ([1999) 
used to estimate stellar temperatures. The star #37345 discussed 
in the text is indicated with an empty triangle. 



was recognized early on as a possible cause of the heavy-mass 
loss underlying the formation of EHB stars, but non-canonical 



mixing phenomena had to be invok ed (von Rudloff et al. 1988; 
Denissenkov & VandenB erg 2003| ISweigart & Mengel 11979 : 
ISweig art 1997). T he discovery of multiple stellar populations in 
GCs (iBedin et alJ 12004 iPiotto et al1l2007l) opened a new fron- 
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tier, because the observational results of Pi otto et aD (120051) ap- 
parently imply that the bluer main-sequence (MS) in to Cen is 
most probably populated by He-enriched stars. The currently 
preferred scenario is that these objects constitute a second stel- 
lar generation that formed fro m material polluted by either 
interm ediate-mass AGB stars (ID' Antona et al. 12001 Renzini 
2008), or rapidly rotating ma ssive MS stars (lMaeder & Mevnet 
2006; iDecressin et alj 120071) . The models of multiple popula- 
tions with different helium abundances successfully reproduce 
both the MS splitti ng and th e multimodal HB morphology of 
both oj Cen (iLee et al.l [2005b and NGC 2808 dD' Antona et alj 
12001 . The He-enhancement thus represents a promising model, 
and alternative to the binary scenario, for the formation of EHB 
stars in GCs. 

In this paper, we present our results of a search for EHB 
binaries in NGC 2808. Preliminary results, pointing to a close 
binary fraction h i gher th an in NGC 6752, were pr e sente d by 
IMoni Bidin et alj d2008bl) and iMoni Bidin & Piottol d2010h . In 
this context, this cluster is a key object because it is notice- 
ably younger than NGC 6752 dDe Angeli et al.ll2005l) . allowing 
the study of the / c -age relation foreseen by the binary scenario. 
Moreover, this cluster represents one of t he greatest successes of 
the He-enhancement scenario, because D' Anto na et a D d2005h 
an dE Orazi et al. ld20ioh were able to model both its multimodal 
HB dBedin etalj|2"00 0) and the multiple MS dPiotto et alJl2007h 
with three stellar populations of increasing primordial helium 
abundance. 

2. Observations and data reduction 

We selected 83 hot HB stars in NGC 2808 from the photo- 
metric catalog of Mom anv et al.l (2003), of magnitude between 
V=17.4, correspond ing approximatively to the Grundahl jump 
dGrun dahl et ali i 19991) . and the limiting magnitude V=19.5 im- 
posed by program feasibility. In Table |4] we give the IDs, 
the coordinate s , and the photometric data of the targets from 
Momanv et al. (2003), and their location in the cluster CMDs 
shown in Figure Q] 

Twenty-four single spectra of target stars were collected 
in fourteen epochs between January 11 and March 24, 2007, 
with the GIRAFFE spectrograph at the VLT-UT2 telescope, in 
both visitor and service mode. The instrument setup H7A pro- 
vided high-resolution (R=18 000) spectra centered on the Hg 
line. The temporal sampling was carefully planned to maxi- 
mize the detection probability of binaries with any period up to 
100 days. Exposures were acquired in pairs to be later summed, 
except when bad weather conditions forced us to stop observa- 
tions after the first frame. We thus collected fourteen epochs of 
data over 2.5 months. The log of the observations is given in 
Table[T] where each exposure is identified with a unique ID, and 
the epoch at the middle of the acquisition period is indicated, 
along with the exposure time and the observing mode (v=visitor, 
s=service). 

Data were reduced with the dedicated CPL-based pipeline 
available at the ESO web site. Because of the extremely low sig- 
nal collected for the hottest targets, we performed many trial re- 
ductions to find the choices and parameter sets that maximized 
the output quality. The frames were de-biased and flat-fielded 
with standard procedures based on the frames collected within 
the standard calibration plan. The dark current was found to be 
non-negligible only along the top edge of the CCD, not used in 
our work, and no dark correction was applied to avoid the corre- 
sponding decrease in S/N by 10-15%. We gave particular atten- 
tion to the wavelength calibration (wlc), whose defects can easily 
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Fig. 2. Normalized spectra of four target stars, obtained shifting 
all the collected spectra to laboratory wavelength and summing 
them. The star ID and temperature is given. The spectra were 
vertically shifted to avoid overlap. 



Table 1. Log of the observations. 



Exposure ID 


Epoch 


U'xp 


Mode 




JD-2450000 


s 




1-1 


54111.26437 


2x 3000 


V 


1-2 


54111.33528 


2x 3000 


V 


2-1 


54112.25032 


2x 3300 


V 


2-2 


54112.32817 


2x 3300 


V 


3-1 


54114.24152 


2x 3000 


V 


3-2 


54114.29434 


lx 3000 


V 


4-1 


54116.25362 


2x 3375 


V 


4-2 


54116.31284 


lx 3375 


V 


5-1 


54127.16030 


lx 1980 


s 


5-2 


54127.25886 


2x 2770 


s 


5-3 


54127.33708 


2x 2770 


s 


6-1 


54182.10276 


2x 2770 


s 


6-2 


54182.16923 


2x 2770 


s 


6-3 


54183.06863 


lx 2770 


s 



affect the radial velocity (RV) measurements. The goodness of 
the wlc was checked by analyzing the spectra of the lamp fibers 
acquired simultaneously with target stars. This reduction step 
was particularly problematic, because we found that running the 
complete wlc routine resulted in an incorrect solution, with a 
deviation from the correct one that increased with wavelength 
and fiber number, up to 10-15 km s _1 . We therefore adopted the 
standard solution for the H7A setup, included in the instrumen- 
tal package downloadable from the GIRAFFE web site, allow- 
ing the pipeline to use the lamp fibers to find rigid shifts and 
changes in the spectral geometry on the chip. After the final ex- 
traction, the lamp fibers showed only small random deviations 
from laboratory wavelengths (0.3 km s _1 rms). This wlc error 
is small compared to uncertainties in the RV measurement, and 
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Fig. 3. Difference between the RV measured in the first frame 
(frame 1 - 1 of Table [TJ and the weighted average for each star, 
plotted as a function of the fiber number. Only the 53 brightest 
stars, with measurement uncertainty lower than 5 km s , are 
plotted. 
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Fig. 4. Mean error of each star as a function of effective temper- 
ature. 



can be safely neglected in the final error budget. Finally, sci- 
ence spectra were extracted using both an optimum algorithm 
(lHornelll986l) and a simple sum. We found that these two meth- 
ods were in general equivalent and the choice did not alter the 
results, but in some noisy spectra one or the other returned more 
precise measurements. This was probably due to small cosmetic 
defects or noise spikes being treated differently by the two algo- 
rithms. We therefore preferred optimum-extracted spectra, but 
we opted for a simple sum in the few cases in which this clearly 
returned smaller RV errors. The background flux was estimated 
by averaging nine fibers allocated to the sky and, after subtract- 
ing their mean spectrum from those of the targets, we checked 
that the weak interstellar emission in the core of the H^g line 
had been effectively removed. The spectra were then trimmed 
to retain only the central region (4780-4930 A), and we nor- 
malized them fitting a linear relation to the continuum on both 
sides of the Hp line. We verified that a higher order polynomial 
was not required in the normalization, as there was no appre- 
ciable change in either the fitted function and or the results. As 
a final step of the reduction, the spectra forming a pair of expo- 
sures (see Table[T} were added. Some example spectra are shown 
in Figure [2] for two stars at the edge of the temperature range 
(12 700 and 21 200 K) and one of intermediate temperature, plus 
the star #37345, discussed later. The presented spectra are the 
sum of all the spectra collected for each star, after shifting them 
to laboratory wavelengths. 

The obse rved HB w a s fitte d with the zero-age HB model 
(ZAHB) of ICassisi et all d 19991) with metallicity [Fe/H]=-1.10 
dCarretta et al.l l2006) to derive a temperature scale along the HB. 
The procedure was not straightforward using the ( U - V) color, 
and uncertainties remained in the determination of the required 
distance modulus and reddening. The se problems could be due 
to the use of the U band, because Dalessand ro et alJ (2010) 
showed that, at shorter wavelengths, it is impossible to fit the 
HB of NGC 2808 with one single population of fixed helium 
content. Therefore, the fit was performed in the V versus (here- 
after vs.) (B - V) plane, where the closest match was found by 
assuming that (m-M)y-15.7 and E(B - V)=0.15, in good agree- 
ment with Bedin et a il (120001) . The fit is shown in the lower panel 
of Figure Q] The temperature of each target was then estimated 
from the point of the model ZAHB closest to the observed po- 
sition. Varying (m-M)v and E(B - V) between values that still 
gave a reasonably good fit, we estimated that the uncertainty in 
the temperature should be on the order of 10%. The temperature 
derived for each star is given in Table |4] 

In the selected sample, there are eight stars hotter than 
20000 K and, following the canonical definition, they can be 
considered EHB stars. However, in the CMD of NGC 2808 
there is no underpopulated region in correspondence to this tem- 



perature as in, for example, NGC 6752 and M80, where the 
EHB i s separated from cooler HB stars. In contrast. ISosin et all 
(119971) identified two clear gaps in NGC 2808 at about V=18.4 
(-16500 K) and V=20 (-25 00 K), called Gl and G2 re- 
spectively, by Bedin et al. (2000). The cluster blue HB is thus 
divided into three sections, called EBT1 (T e ff <16500 K), 
EBT2 (16000<T eff <25 000 K), and EBT3 (T eff >25 000 K). 
Following this scheme, the sample contains 62 EBT1 and 21 
EBT2 stars, but no star in the faint EBT3 group was observed. 

We note that the temperature associated with the gaps was 
obtained by fitting a canonical mo del to the observed H B, and 
the spectroscopic measurements of M oehler et al.l (120041) agree 
with this temperature scale. However, Dalessandro et al. (2010) 
proposed a new scale, based on a multi-population model of 
three stellar generations with increasing helium content. In their 
calculations, the Gl and G2 gaps approximatively coincide with 
the canonical start (20000 K) and end (31 000 K) of the EHB, 
which is thus separated from cooler stars by an underpopulated 
region even in the CMD of this cluster. 

3. RV measurements 

The RVs were measured by means of the cross-correlation 
(CC) technique dTonrv & Davis! 1 19791) implemented in the fx- 
cor IRAlfl task. The CC was restricted to the spectral interval 
4840-4880 A, i.e. the Hg line with its full wings. The templat e 
was extracted from the synthetic library of Mun ari et al.l d2005). 
Experiments with model spectra of different temperature, grav- 
ity, and metallicity showed that a change in these parameters 
does not affect the results, while the errors are very sensitive to 
this choice, and the templates with a narrower Hp usually provide 
smaller errors for all the stars. T hese conclusions ar e also sup- 
ported by the detailed analysis of Mors e et al 1 (119911) . Variations 
in the shape of the template line are indeed not expected to shift 
the center of the CC function (CCF), but can affect the resulting 
uncertainties. We also found that the inclusion of weak metallic 
lines in the synthetic spectra enhanced the errors without adding 
real information, because they were either absent or not visible 
in the noisy target spectra. Consequently, the template finally 
adopted for all the stars was a synthetic profile of the Hp line, 
obtained by removing all weaker lines from a synthetic spec- 
trum at 20000 K, log(g)=5, and cluster metallicity. It must be 
noted that the RV measurements are unaffected by any bias or 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 5. Upper panel: histogram of the distribution of the residu- 
als with respect to the weighted average of each star, in unit of 
the associated error (Equation [T}. Overplotted to the histogram, 
a Gaussian with cr-0.96 and centered in 0.01 is shown. Lower 
panel: probability plot of the same residuals of the upper panel. 
The straight line, with intercept 0.01 and slope 0.96, indicates 
the least squares fit of the data. 



uncertainty in the estimate of the target temperature, because the 
choice of the template was independent of this temperature and, 
as discussed, this choice did not alter the results. 

The RV was determined by fitting the core of the CCF with a 
Gaussian profile. While the procedure was straightforward for 
cool stars, it became problematic for the fainter ones, whose 
spectra were much noisier down to signal-to-noise ratio (S/N)f»3 
in the worst cases. Measuring RVs at such low S/N with only one 
wide line is a challenging task, and we with experimented dif- 
ferent methods to improve the accuracy of the measurements : 
the application of a narrow Fourier filter rtBrault& Whitelll97 lh 
sometimes helped, while the S/N of very noisy spectra was of- 
ten increased by degrading the resolution by a factor of 2-3 and 
then rebinning accordingly. Although a lower resolution affects 
the precision that can be achieved, the noise is by far the domi- 
nant source of uncertainty when the stellar flux is very weak. In 
some cases, we used only one spectrum of a pair, because there 
was a great difference in the spectral quality of the two, and the 
addition of the noisier one degraded the resulting measurement. 
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BO 90 100 110 
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Fig. 6. Distribution of the absolute RVs of the observed 
stars. The curve shows a Gaussian centered on the mean 
value (93.2 km s _1 ) and cr equal to the observed disp ersion 

5.7km s" 1 ). The dotted line indicates the cluster RV from Hams 

1996T) . 



We did not find a unique scheme for returning the optimal re- 
sults for all the spectra, but we employed in each case the pro- 
cedures (filtering, rebinning, use of a single spectrum) providing 
the best results in terms of shape and noise of the CCF, height 
of its central peak, and goodness of its Gaussian fit. These pro- 
cedures were designed only to reduce the errors, and the results 
were stable when experimenting with different combinations of 
procedures and involved parameters, varying by no more than a 
few km s _1 . When this was not the case, the measurement was 
judged unreliable and excluded, as when the results were sensi- 
tive to changes in either continuum normalization or extraction 
algorithm. 

3.1. Systematic errors 

After the correction to heliocentric RVs, we verified that no 
systematic error was present in the results. First, we checked 
the zero-point of each frame by averaging the RVs of the 53 
brightest stars, excluding measurements with errors larger than 
5 km s~ 1 . We thus derived the corrections to reduce each frame to 
the same zero-point, although they were lower than 1.5 km s _1 , 
i.e. well within the typical errorof the 53 stars (3.5-4 km s~'). We 
then plotted, for each frame, the residual of each star with respect 
to its weighted-averaged RV as a function of the fiber number, to 
check for the presence of a sys tematic effect that varied with po- 
sition on the CCD, as done by Moni Bidin et al. ( 2006, Figure 5 
and 7). The plot relative to the first frame is shown in Figure [3] 
as an example. The average value of these residuals was always 
lower than 0.1 km s _1 , indicating that any offset between expo- 
sures was correctly removed in the previous step. Moreover, the 
linear and third-order fit never differed from zero by more than 
1 km s _1 , proving that there is no residual trend in the measured 
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Fig. 7. Upper panel: Probability of a false detection for the "rel- 
ative" criteria (empty dots and dotted line) and the "absoute" 
criteria (full dots and line) defined in the text, assuming a 3<x 
threshold for both. Lower panel: False detection probability for 
the same criteria, as a function of the threshold, for ten measure- 
ments. 



RVs, and the results are free from systematics well beyond the 
typical random errors of 3-4 km s _1 . 



3.2. Errors 

Any search for binaries by means of multi-epoch RV measure- 
ment is based on the comparison between the observed varia- 
tions and the uncertainties. A precise definition of the errors is 
therefore of fundamental importance for the correct interpreta- 
tion of the results. 

The uncertainty associated with each measurement was de- 
fined as th e quadratic sum of th e error of the CC technique (as 
defined bv lTonrv & Davisll979l) . the wlc error defined in Sect. [2] 
and the uncertainty introduced when correcting the zero-point 
of each frame (Sect. 13. 11 1. This last quantity was defined as the 
rms of the corrections applied to each frame, i.e. the scatter in 
the zero-point of the frames around the mean value. The CC 
error completely dominates the error budget, being typically 2- 
3 km s _1 for bright stars and increasing with temperature up to 
5-8 km s _1 for the faintest targets. Thus, the quadratic sum of 
the wlc and zero-point errors (0.3 and 1.1 km s _1 , respectively) 
makes a negligible difference for all but the brightest targets. The 
mean error for each star is plotted in Figure 0as a function of ef- 
fective temperature. The errors constantly increase with the tem- 
perature along the HB, reflecting the decreasing S/N of the col- 
lected spectra. The mean value of EBT1 targets (T e jf < 17 000 K) 
ranges from 2.5 to ~4.5 km s _1 , with the exception of a couple 
of deviant stars, with an average value of 3.5 km s . Hotter stars 
show a higher scatter and a steeper gradient, with a mean value 
of 6.0 km s _1 . 



To test the reliability of the errors, we analyzed the distri- 
bution of the residuals (A,) of the i-th measurement (RV,) with 
respect to the weighted average (RV), in units of the associated 
error cr RV ,i 



RVi - RV 

0"RV,i 



(1) 



If the observational uncertainties are the only cause of the vari- 
ations, and the errors (Trv^ represent them well, the residuals 
would follow a normal distribution centered on zero and with 
unit dispersion. This analysis can easily be spoiled by RV vari- 
able stars, which add large values to the wings of the distribu- 
tion. We therefore excluded all the stars showing one or more 
residuals A,- > 3.5, as we later demonstrate that no variation 
above this threshold should be expected from random errors 
alone (Sect. 14. il l. The resulting distribution of A, is shown in 
the upper panel of Figure [5] Both the mean value (0.01) and 
the standard deviation (0.96) confirm that the errors are indeed 
well defined, as they account exactly for th e random variations. 
The probability plot (iLutz & Hanson! 19 92) of these data reveals 
that they follow a normal distribution, as they are aligned along 
a straight line of slope 0.96 and intercept 0.01. This would have 
not happened if, for example, the exclusion of stars with A, > 3.5 
had produced a too narrow cut of the wings of the distribution. In 
conclusion, the analysis reveals that the errors are reliable, rep- 
resenting well the random uncertainties, and that they can safely 
be used in a statistical analysis. 

3.3. Absolute RVs and cluster membership 

The absolute RVs were calculated by weight-averaging all the 
measurements for each star, and the results are given in col- 
umn 7 of Table [4] The histogram of the distribution is shown 
in Figure [6] The m ean value is R y ar ,s=93.2 km s _1 , in excel- 
lent agreement with Harris (1996, February 2003 Web version, 
93.6 km s~'), and the observed dispersion is 5.7 km s -1 . In 
Figure [6] we also overplot a Gaussian curve centered on R V a b s 
with cr-5.1 km s _1 . The histogram is well described by a normal 
distribution, with no deviating points, hence all the targets can 
be considered bona fide cluster members. The only exception is 
the star #37345, whose spectrum is shown in Figure [2] This ob- 
ject was immediately recognized as peculiar, because of its very 
broad Hg line: rapidly r otating HB stars hotter than ~1 1 500 K 
have never been found dRecio-Blanco et al.1 12001 lBehril2003l) . 
while they are very common among hot MS stars. We note that 
this object is much redder than the other HB stars in terms of 
U — V color (of Figure [TJ, but merges with the HB population in 
the V-(B - V) diagram (lower panel of the same Figure). The RV 
measurements was highly uncertain because, as a consequence 
of the line broadening, the fit of the wide peak in the CCF was 
problematic. When it converged, we derived RV*150 km s , 
very different from the cluster value. We therefore conclude that 
the star #37345 is most probably a background field MS B star. 

4. Analysis 

Eleven stars were not considered in our analysis, because their 
spectra were either too noisy for reliable measurements or 
strongly contaminated by the emission lines of nearby lamp 
fibers. All these targets have no data in the last two columns 
of Table [4] The majority of targets were not allocated the same 
fiber in both GIRAFFE plates, as these are not equivalent be- 
cause of broken fibers. As a consequence, the spectra of six stars 
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were free of lamp contamination when observed with one plate, 
while their spectra collected with the other plate were damaged. 
Moreover, star #45980 had no fiber allocated in one plate. These 
seven targets were studied as the others, but were excluded from 
our statistical analysis because, with only half of the spectra 
with respect to the other stars, their temporal sampling is very 
different. Finally, one target was found to be a foreground star 
(Sect. l3.3t . In conclusion, with nineteen stars lost or excluded, 
the resulting sample comprises 50 EBT1 and 14 EBT2 stars, out 
of which six EHB targets have T eff >20 000 K. 

The measurements for the four spectra that were not the sum 
of two exposures, i.e. the frames 3-2, 4-2, 5-1, and 6-3 (TableQ]), 
were in general less reliable, and they were often excluded for 
faint stars. In the upper panel of Figure [8] we plot the decrease 
in the binary detection efficiency, as a function of period, if these 
four frames are not considered. Their exclusion clearly does not 
affect the survey much, because the probability of detecting a 
true binary in general decreases by less than 3%, except for two 
sensitive periodicities (1 day and ~18 days) where the loss is 
about 10-15%. Therefore, we finally excluded them from the 
statistical analysis. The measurements for these lower-quality 
frames were indeed always more uncertain, and were often ex- 
cluded anyway, affecting the uniformity of the measurements. 
We then consider only the remaining ten epochs for each star. 

4.1. Detection threshold and detection probability 

The identification of binary candidates requires criteria to de- 
fine when the observed variations can be considered significant. 
The criteria must satisfy two desiderata: the probability of a false 
detection (Pfaise) must be negligible, while the probability of de- 
tecting a true binary (Pdet) must be as high as possible, within 
the limitations imposed by the temporal sampling and the obser- 
vational errors. Usually the first point is satisfied if the statistical 
expectation is less than one false detection in the whole survey, 
which in our samples of 64 targets means Pf a i se <0.014. A com- 
promise between these requirements is often needed, because 
more stringent criteria reducing Pf a i se also reduce the efficiency 
of detecting genuine RV variables. 

The parameter Pd e t was estimated as a function of orbital pe- 
riod, generating 2500 synthetic binaries in circular orbits of pe- 
riod p, comprising two stars of 0.5 M Q , uniformly distributed in 
the sin(/)-6> space, where i is the angle of inclination of the or- 
bital plane and 6 is the orbital phase. These systems were then 
"observed" with the same temporal sampling of our survey, and 
each star satisfying the criteria under analysis represented a de- 
tection. The fraction of detections over the whole sample thus 
indicated the efficiency of the survey for systems of period p. 
Pfaise was calculated by simulating 100000 sets of N measure- 
ments drawn from a normal distribution centered on zero and of 
unity dispersion. Each set represented the observations of a star 
with constant RV, affected by only random errors normalized to 
unity, and N-10 for our survey. A false detection was claimed 
for each set of measurements satisfying the criteria for binary 
detection, and Pf a ] se was estimated as the fraction of false detec- 
tions over the whole sample of 100000 attempts. 

In previous investigations of HB stars in GCs, the criteria 
adopted for the detection of a binary candidate was that the mea- 
sured RV variation was larger than 3<x, where <x was either the 
error i n the s pectral shift among two epochs dMoni Bidin et al.l 
2006, 2009a), or the quadratic sum of the errors in RV mea- 




surements (Moni Bidin et al. 2008a). This strategy worked well 
for surveys based on 4-5 epochs, but is unsuitable for our work: 
the upper panel of Figure [7] shows that the probability of a false 



1 10 

Period [days] 



Fig. 8. Upper panel: difference of the binary detection proba- 
bility, as a function of period, when including or not the four 
spectra which are not the sum of two exposures (frames 3-2, 4- 
2, 5-1, and 6-3, see Table [TJ. For the calculation we assumed 
an "absolute" criteria with threshold 3.2cr, the same adopted in 
the present work, and a mean observational error of 3.5 km s _1 . 
Lower panel: detection probability as a function of period, calcu- 
lated as described in the text. The full and dotted lines indicate 
the results for the "absolute" criteria with 3.2<x threshold and 
the "relative" one with 3.8cr criteria, respectively, assuming the 
typical mean error (3.5 km s" 1 ) of stars cooler than the Gl gap 
at ~17 000 K. The results for the "absolute" criteria assuming a 
mean error of 6.0 km s -1 , characteristic of hotter stars, are shown 
with a dashed line. 



detection rapidly increases with the number of measurements, 
and in our 10-epochs survey we should expect 8% of non-binary 
targets (~6 stars) to violate the threshold due to random errors 
only. The method could be applied even here, but to reduce the 
expected false detections to less than one star we should increase 
the threshold up to 3.8cr, as indicated by the lower panel of 
Figure [7] i.e. ~25-30 km s _1 for the hotter stars. In the same 
Figure, we also show Pf a i se for an alternative criteria, which we 
define as "absolute" (while the older one is called "relative"): a 
star is flagged as binary candidate if one (or more) measurement 
differs by more than a certain threshold, in units of its associ- 
ated error, from the weighted average. The plot reveals that the 
reliability of this criteria is poorly affected by the higher num- 
ber of measurements, and that in our survey the threshold must 
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Fig. 9. Maximum residual in units of the error (Equation [T] for 
each star, as a function of temperature. The horizontal line in- 
dicates the adopted threshold for binary detection. Empty dots 
represents stars with only partial temporal coverage, that are ex- 
cluded from the statistical analysis. The vertical dotted line indi- 
cates the approximate position of the Gl gap. 



be fixed to 3.2cr to have P&ise <0.7%, i.e. less than 0.5 false de- 
tections expected among our 69 targets. In the present work, we 
adopt both this criteria and this threshold for the detection of bi- 
naries. However, the results are completely equivalent to those 
for a "relative" criteria with threshold of 3.8cr, because Pf a i se is 
the same, and Pd e t is very similar to within a few percent at all 
periods (see lower panel of Figure [8). 

The efficiency of the survey resulting from the adopted cri- 
teria, i.e. the probability of detection as a function of the orbital 
period, is shown if Figure [8] We considered two possibilities: a 
mean error of 3.5 km s _1 , typical of the targets brighter than the 
Gl gap (T eff <17 000 K), and 6.0 km s" 1 for the hotter EBT2 
stars. In the first case, the survey can detect binaries with pe- 
riods up to 400 days and, although the probability of detection 
drops below 50% for periods longer than 50 days, still more than 
one-third of the binaries with period 50-200 days can be found 
with our observations. For hotter stars, the larger errors decrease 
at all periods the efficiency, which becomes negligible (<20%) 
for P>50 days. In any case, the presence of HB binaries with pe- 
riods longer than 10 days can be tested for the first time in a GC, 
and these wide systems are poorly studied even among field sdB 
stars. 

4.2. Detected binaries 

The results of the search for RV variable stars are summarized 
in Figure [9] where for each star we plot the maximum residual 
Amax defined in EquationQ] This value is also given in column 8 
of Table |4] Five targets show one or more significant deviations 
from the weighted mean, and they must be considered candidate 
binaries. The cooler one (star #32670) is one of the seven targets 
with only partial temporal coverage mentioned at the beginning 
of this section, and is excluded from the statistical analysis. The 
other four binary candidates are located in the hotter half of the 
temperature range, and are all EBT2 stars fainter than the Gl 
gap. All but one of them (#9715) are hotter than 20000 K, thus 
fall within the canonical definition of EHB stars. Their measured 
RVs are shown in Figure [10] 

We note that, although the hotter stars were excluded from 
the error analysis of Sect. 13. 21 the stars with no variation above 
the threshold show no gradient in A max with temperature, indi- 
cating that the errors are well-defined across the whole range. 
Moreover, there is a gap in the distribution of A max between these 
stars and the five binary candidates, which reinforces the idea 
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Fig. 10. RV measurements of the four binary candidates. The 
dotted line indicates the weighted mean. 



that they are outliers. This confirms that the adopted threshold is 
indeed a good separation between statistical random errors and 
real RV variations. 

4.3. Period estimate 

Our RV measurements are too few in number, and affected by 
too large errors, to attempt any estimate of the orbital periods of 
the binary candidates. However, our survey can detect systems 
with periods up to 400 days (Sect. l4.ll i. while previous inves- 
tigations of GCs and the most extensive surveys among field 
stars were limited to close binaries with periods shorter than 
10 days. This implies that our results are not directly compa- 
rable to those works without a minimum classification of the 
discovered systems. Hereafter, we define close binaries as sys- 
tems with p <10 days. They have been the common targets of 
previous investigations, and should have undergone one or two 
CE phases. We also define intermediate-period binaries as ob- 
jects with 10< p <400 days. We thus leave the definition of wide 
binaries to those with pe riods longer than 400 days, which is the 
typical product of RLOF dHan et alj|2002l) . a kind of system un- 
detectable in our survey. 

To help classify the binary candidates, we can rely on the 
adopted temporal sampling: the observations of the first four 
closely-spaced nights were designed to detect short-period sys- 
tems, while the later fifth and sixth epochs were planned to help 
us discover longer periodicities. If our observations had been re- 
stricted to the first two nights only, the stars #9715 and #55759 
would still have been detected. As can be seen in the upper panel 
of Figure QT| systems with periods longer than 5 days would 
hardly have been found in such a two-night survey, and all bina- 
ries with p >10 days would have passed completely undetected. 
The period of these two binary stars must therefore be shorter 
than 10 days, and most probably even shorter than 5 days. These 
objects can thus safely be classified as close binaries. Figure [TT] 
also reveals that ~80% of the close binaries should have been 
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Period [days] 

Fig. 11. Upper panel: detection probability as a function of pe- 
riod, calculated as described in the text, consideting only the first 
four (full line) or two (dotted line) nights. Lower panel: proba- 
bility, as a function of period, that an EHB binary would show a 
significant RV variation only in the fourth night of our observa- 
tions, as found for the star #7700. 



detected during the first four nights. The remaining 20% mainly 
consists of nearly face-on systems showing maximum RV varia- 
tion near or below the detection threshold, as indicated by the in- 
clusion of the other two epochs not noticeably increasing the sur- 
vey efficiency (compare Figs. [TTIand[8l>. The star #9519 shows 
variations above the threshold, but not in the first four nights, and 
therefore can safely be classified as an intermediate-period bi- 
nary of period longer than 10 days. The classification of the last 
binary candidate (star #7700) is less straightforward, because it 
passed undetected during the first three nights but it had a signif- 
icantly different RV during the fourth one. To clarify its nature, 
we calculated the probability that a system of a given period be- 
have as observed, i.e. the fraction of the synthetic binaries de- 
fined in Sect. l4.1l that would show a significant variation in the 
fourth epoch only. The results are plotted in the lower panel of 
Figure QT| The period is most likely between 5 and 20 days, 
although shorter values are not impossible, and the most prob- 
able period is ~1 5 days. Taking into account that, as argued by 
M oni Bidin et alj d2008al) . restricting the definition of close bi- 
naries to p <5 days would not make a great difference, we can 
classify this object as an intermediate-period binary. 



of N targets were detected. The calculation requires informa- 
tion about the underlying period distribution, which is unknown. 
A Gaussian distribution in log (p) centere d on zero has been 
proposed for close field sdB systems (e.g. iMaxted et alJ 12001: 
Napiwotz ki et al.1 1200 4). while wider systems are still poorly 
studied. We assume a fl at distribution at all periods, because 
M oni Bidin et al . (2008a) demonstrated that the results do not 
change by more than a few percent if other reasonable, but not 
proven, assumptions are made. 

The binary fraction will be estimated separately for the 
EBT1 and EBT2 sections of the HB, which have different sam- 
plings (Figure [TJ and detection probabilities (Figure |8}. These 
two groups of stars are separated by an underpopulated region 
in the CMP, and cou ld represent different stellar populations 
(iD'Antona et al.ll2005l) . Stars hotter than 20000 K will be con- 
sidered as EHB stars, a sub-group of the EBT2, although this 
could be only a formal distinction since no known discontinu- 
ity separates them from the other EBT2 stars in the CMD. It 
must also be recalled that this definition strongly depends on the 
model fixing the tem perature scale along the H B, and the multi- 
population model of iDalessandro et ail d2010l) EBT2 and EHB 
actually coinciding. We also distinguish between the fraction of 
close and intermediate-period binaries (f c and /; p , respectively). 
The analysis of fi p is restricted to p— 10-200 days in the EBT1 
and to p=10-50 days in the EBT2 and EHB, because the effi- 
ciency of the survey rapidly decays for longer periods. In the 
first case, this choice does not affect the results, because no bi- 
nary was detected and therefore the solution Nb-0 can be ap- 
plied to any period interval of interest. In contrast, this restriction 
could lead to wrong results for the EBT2/EHB because, if one of 
the two candidate intermediate-period binaries had p >50 days, 
/ip would be overestimated. However, we have already demon- 
strated that the star #7700 has p <30 days (see Figure fTTT i. and 
a similar analysis for the star #95 19 suggests that the probability 
that p >50 days is negligible even for this star. 

4.4.1. EBT1 

In the EBT1, we have N-50, and Nb=0 for both close and 
intermediate-period systems, with J=0.86 and 0.48, respectively. 
The resulting probability curves for the two binary fractions are 
shown in the upper panel of Figure Q~2] The curves are narrow- 
peaked because the efficiency of the survey is high and the sam- 
ple is sufficiently large, with the curve for /; p being wider be- 
cause of a lower detection probability. As we failed to detect bi- 
naries in the EBT1, the most probable value is always zero. The 
probability of our observations drops below 10% for f c >5% and 
/ip >9%, and below 5% for f c >6% and / ip >1 1%. These can be 
adopted as upper limits at the 90% and 95% confidence levels, 
respectively. 



4.4. Binary fraction 

The probability of detecting Nb binaries in a sample of N stars 
is given by the expression 



N! 



(N-N B )!N B : 



(d/) NB (l-d/) 



N— Nb 



(2) 



where / is the underlying binary fraction and d is the aver- 
age probability of detection weighted by the period distribu- 
tion. This equation can also be interpreted as the probability 
that the (unknown) binary fraction is /, when Nb systems out 



4.4.2. EBT2 

In the EBT2, we have A^=14, and Nb=2 for both close and 
intermediate-period systems, with J=0.74 and 0.43, respectively. 
The curves of probability for the two binary fractions are shown 
in the middle panel of Figure [12] The most probable estimate of 
f c is 19%, and the detection probability is high enough to obtain 
a relatively narrow curve: the 90% and 95% confidence level in- 
tervals are 6-42% and 4-49%, respectively. Stronger constraints 
could have been obtained only with a larger sample, but this is 
a typical observational limitation for GCs, because hot, isolated 
stars outside the inner crowded regions are few in number. The 
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Fig. 12. Curve of probability for the close binary fraction f c 
(dotted lines), and the intermediate-period binary fraction /; p 
(full lines). Upper panel: EBT1 stars; middle panel: EBT2 stars; 
Lower panel: EHB stars. 



best estimate for /j p is 33%, but the large errors affected the ef- 
ficiency at detecting larger periods, and the probability curve is 
very wide. The 90% and 95% confidence level intervals are 11- 
72% and 7-85%, respectively. These results can serve as lower 
limits, since very low intermediate-period binary fractions can 
be excluded. 



4.4.3. EHB 

In the EHB, we have N-6, N B =h and d-Q.lA for close bina- 
ries, while Nb=2 and J=0.43 for wider systems. The sample 
is too small to allow us to place strong constraints on the bi- 
nary fractions, in particular for /j p where the detection probabil- 
ity is low. Therefore, we analyze this sub-group of EBT2 stars 
only for completeness. The more probable values are / c =23% 
and /; p =77%, and the close binary fraction is contained within 
the interval 3-66% at the 90% confidence level, or 2-78% at the 
95% level. Only lower limits can be obtained for fi p , because the 
probability curve does not decrease much after its maximum: at 
the usual 90% and 95% confidence levels, f ip >24% and 16%, 
respectively. 

5. Discussion 

The results for the close and intermediate-period binary fraction 
(/ c and fip) obtained in Sect. l4.4.Tll4.4.3l are resumed in Table [2] 
where we present the most probable estimate (/ max X an d the 90% 
and 95% confidence level intervals (P(/)>10% and P(/)>5%) 
separately for the EBT1, EBT2, and EHB stars. Despite our ef- 
forts to obtain the most precise measurements, the very low spec- 
tral quality of the fainter stars strongly affected our results: the 
detection efficiency was reduced by the large errors, especially 
for the longer periods, and one-fourth of the observed EBT2 stars 
had to be excluded. As a consequence, the binary fractions are 
not very well determined, as reflected by the large confidence 
intervals for the EBT2 and EHB groups. However, the statistical 
analysis still permits us to draw many interesting conclusions. 

The calculations of Sect. 14. 41 are based on the fraction of 
detected binaries relative to the number of stars observed in 
each group. Some systematic error in the estimated tempera- 
tures could thus alter the results, in particular in the presence 
of stars evolving off the EHB toward higher luminosities, whose 
temperature could be underestimated. The temperature and lu- 
minosity of the hotter stars are degenerate in the V vs (B - V) 
diagram, because the bluer section of the HB is almost vertical. 
Small luminosity variations or temperature errors should not af- 
fect the definition of the EBT1 and EBT2 samples, because it 
is based on the presence of clear photometric gaps and not on 
temperature boundaries. Stars ascending the asymptotic giant 
branch can be identified in the V vs (U - V) plane, where the 
temperature-luminosity degeneracy is much less pronounced, 
and EBT3 stars, evolving directly to the WD cooling sequence, 
should be at all stages fainter than our magnitude limit. As a con- 
sequence, targets assigned to the EBT2 sample should be cor- 
rectly classified, while two EBT1 targets brighter and/or bluer 
than the main HB population (see he upper panel of Figure [TJ 
could be evolved post-HB stars originating from the EBT2 pop- 
ulation. One of them was excluded from the statistical analysis 
for other reasons (Sect.|4), while the other one shows no sign of 
binarity. Assigning it to the EBT2 group, we would have N-15 
in Sect. l4.4.2l causing a change of about 2% in the values given 
in Table|2] The results for the EBT1 group are even less sensitive 
to the change, because of the larger sample. In conclusion, the 
uncertainties in the definition of the EBT1 and EBT2 samples 
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have a negligible impact on the results. In contrast, even small 
systematics in the temperature scale can strongly affect the sta- 
tistical analysis of the EHB group, defined only by a temperature 
limit, because of the small number of observed targets. For ex- 
ample, when N is varied by +1, / C;lnax and /j Pima x can vary by 
+5% and +15%, respectively. These changes are not large when 
compared to the uncertainties given in Table [2] and the general 
conclusions are unaffected, but this result warns us that the re- 
sults for the EHB group should be regarded as only indicative, 
as already stated in Sect. 14.4.31 

Binaries with periods longer than 10 days have never been 
studied in GCs, and have be en the target of very few su rveys 
even among field stars (e.g. iMorales-Rueda et a l. 2006). The 
most probable estimate of the intermediate-period binary frac- 
tion among EHB stars is very high, but is most probably affected 
by small number statistics as can be deduced from the very wide 
probability curve in the lower panel of Figure [12] However, an 
important result of our investigation is a relatively high (15-20%) 
lower limit for /; p among EHB stars. The estimated close binary 
fraction in all the clusters studied so far is either comparable to 
(NGC 2808, NGC 5986, M80) or much lower than (NGC 6752) 
this value. The value of /i p is most probably higher than / c even 
in the EBT2, because their probability curves are similar but 
p(/i p ) is shifted toward higher values. The most probable val- 
ues indicate that only one-fifth of the EBT2 stars can reside in 
close systems, but that even half of them could be binary systems 
with p <50 days. Although we cannot provide a more reliable 
estimate, this result suggests that binaries wider than those in- 
vestigated so far could play an important role and should deserve 
more attention, particularly where close systems are lacking. 

A second prominent result clearly evident from Table |2]is 
that the 90% confidence intervals in the EBT1 and EBT2 do not 
overlap, for both / c and /; p . This means that the probability that 
the binary fraction has the same value in these two sections of the 
HB is negligible (of the order of 1%): the Gl gap thus appears 
as a discontinuity in the binary fraction, there being a very small 
quantity of binaries among the stars cooler than the gap, and 
about 20-30% for hotter objects. Unfortunately, our data cannot 
exclude / c and /j p monotonically increasing with the temperature 
rather than showing a real discontinuity at Gl. In this regard, we 
note that there is reasonable agreement between the results in 
the EBT2 and in its (hotter) sub-group the EHB, in particular 
for the close binary fraction (compare also the shapes of the dot- 
ted curves in the middle and lower panels of Figure HZt . Thus, 
/ c could be quite homogeneous within the EBT2 population, in- 
stead of increasing with temperature, but no firm conclusion can 
be drawn because of the too wide confidence intervals. Future 
observations should help us to clarify this issue: a sudden in- 
crease in the binarity in correspondence with Gl would strongly 
relate it to the formation of all EBT2 stars, while a / c slowly in- 
creasing with temperature may indicate that the progeny of close 
systems are preferentially hotter. 



In the context of the He-enrichment scenario for the for- 
mation of hot HB stars in GCs, the difference in binary frac- 
tion a mong EBT1 and EBT2 stars contrasts with the investiga- 
tion of ID'Orazi et alj (1201 Ol) in M4, who found a much lower 
quantity of binaries among red giants that displayed evidence of 
chemical enrichment, than among normal ones. They argue that 
this difference is naturally explained by assuming that the sec- 
ond stellar generation formed in a denser environment, where 
more frequent dynamical interactions enhanced the disruption 
rate of binary systems. However, in the multi-population models 
of NGC 2808, the EBT2 is interpreted as the progeny of the lat- 
est and He-richest of the three populations obse rved in the MS 
dD'Antona et all 120051; iDalessandro et al.ll20Tob and, following 
ID'Orazi et al.l (120101) . we would expect the EBT2 to be depleted 
in binaries, at variance with what is observed. If their results 
were confirmed as a general behavior of the chemically polluted 
stars in GCs, our results would argue against the link between 
the EBT2 and the He-enriched stars. However, an alternative in- 
terpretation would be that both the He-enrichment and the bi- 
nary scenarios co-exist in the cluster, as different channels for 
the formation of blue HB stars. In this case, in the EBT2 both 
the progeny of He-enriched stars and products of binary interac- 
tions would be found. This would cause a higher frequency of 
EBT2 stars in the HB with respect to the fraction of He-rich MS 
stars, but binaries in GCs usually represent a minor fraction of 
the entire population, and this difference could pass unnoticed. 
No splitting of the red- or sub- giant branch has been detected 
so far in NGC 2808, but this clu ster shows a strong Na-O an- 
ticorrelation (ICarretta et al.l 12006). which has often been inter- 
preted as a consequence of cluster self-pollution. It would there- 
fore be really instructive to repeat the D'Orazi et al.'s investi- 
ga tion among the three different groups of red giants identified 
bv lCarretta et al.1 d2006) on the basis of their Na/O value, for a 
direct comparison with the values presented in Table [2] 

Previous surveys focused mainly on the close binary frac- 
tion among EHB stars, defined by the temperature boundary of 
20 000 K. This is more than a conventional definition because, 
as noted in Sect.Q] the post-HB evolution of EHB stars is very 
distinct from that of cooler stars. However, the formation mech- 
anism does not necessarily differ, and population synthesis mod- 
els show that sdB stars as cool as 15 000-16 000 K can be formed 
with the same binary channels used to model the field EHB 
star population dHan et al.l2003l) . In addition, massive (0.75 M ) 
EHB stars move to cooler tempera tures in the first stages of the 
post-HB evolution (lHan et alJ2002l) . The previous investigations 
were not limited to stars hotte r than 20 000 K, and the o nly close 
system found in M80, which Mon i Bidin et alJ d2009al) included 
in a statistical analysis similar to ours, is actually cooler than this 
temperature (~18 000 K). For all these reasons, we can safely 
use our estimate of the close binary fraction in the EBT2 popu- 
lation for comparison with both the results for other clusters, and 
the expectations of the binary model for sdB star formation. We 
also recall that the adopted canonical temperature scale might 
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n ot be appropriate, as the E BT2 and EHB coincide in the model 
of lDalessandro et alJ (l2010h . 

The most reliable estimate of f c is low compared to field val- 
ues, as observed in all the GCs observed so far. Even the low- 
est fractions proposed in the literature for the field sdB pop- 
ulation are very unlikely ( ~45%, e.g. iNapiwotzki et aD [2004: 
iMorales-Rueda etal . 2006), and the h ighest ones can safely be 
excluded (~70%. iMaxted et alj 120011) . Howe ver, even the ex- 
tremel y low f c measured in NGC 6752 (4%, IMoni Bidin et all 
2008a) can be excluded with a 95% confidence level. We 
thus find that NGC 28 08 is similar to NGC 5986, where 
M oni Bidin et al. (2009a) estimated that / c =25%, although with 
great uncertainties. The similarity can be drawn even further, be- 
cause both these cluste rs are 2-3 Gyr younge r than NGC 6752 
(|De Angeli et alJl2005l) . While the results of IMoni Bidin etalJ 
d2009al) had to be confirmed, the probability that the close bi- 
nary fraction is equal or very similar in all these three clusters is 
negligible, and the trend of f c with age depicted in Figure 1 of 
M oni Bidin & Piottd(l2010h is likely r eal, following the close bi- 
nary fraction-age relat i on pro posed by [Moni Bid in et al.l d2008a) 
and modeled by iHanl (|2008) in the context of the binary sce- 
nario. On the other hand, NGC 280 8 also confirms the prob - 
lems of this scenario p ointed out by IMoni Bidin et alj ((2009a): 
no model proposed by IHanl (120081) can predict both the very 
low f c of NGC 6752 and a steep increase with decreasing age 
up to f c k2Q% for clusters 2-3 Gyr younger. As can be seen 
from their Figure 3, the proposed solutions do not have a suffi- 
ciently steep gradient, and all the curves with / c (13 Gyr)<5% are 
nearly fiat for old populations, so that / c <10% fo r any cluster 
older than 8 Gyr. Nevertheless, the simulations of IHanl ([2008) 
depend strongly on a set of poorly known parameters. In par- 
ticular, Marsh et al.l (1 19951) deduced a very high CE ejection ef- 
ficiency (acE ~1) from their observations of double-degenerate 
close binaries, b ut a much lower value (0.2-0.3) was proposed by 
IZorotovic et al.l (120101) as a result o f their study of a large sam- 
ple of post-CE binaries. Han (2008) assumed a acE=0. 5-1.0, but 
from their Figure 3 it can be seen that lower values tend to return 
a steeper slope of the / c -age relation, although still predicting too 
high values at 13 Gyr. It is therefore possible that, with a differ- 
ent set of parameters, the binary scenario can be reconciled with 
the observations. 

6. Conclusions 

We have analyzed the radial velocity variations of a sample of 50 
EBTl stars (T eff =12000-17000 K) and 14 EBT2 (T eff =17000- 
22 000 K) stars in the horizontal branch of NGC 2808, which 
were spectroscopically observed in ten epochs spanning a tem- 
poral interval of ~75 days. We detected two close binaries (pe- 
riod p <10 days), and two wider systems (p >10 days and 
p=5-2Q days, respectively). All these object belong to the EBT2 
population, while no binary was detected among the cooler 
EBTl stars. We estimated the most probable fraction of close 
(p <10days) and intermediate-period (p=10-200 days) binaries, 
f c and fp, plus the ranges corresponding to a 90% and 95% con- 
fidence level, respectively, in both the EBTl and EBT2 groups 
of stars. Although for hot stars both the sample and the survey 
efficiency were reduced by the low S/N of the spectra, we were 
able to draw some important conclusions: 

- We found a relatively high lower limit (15-20%) for f v 
among EHB stars (T e jf >20000 K). The intermediate-period 
binary fraction is most probably higher than f c even for 
EBT2 stars, where only one-fifth of stars could reside in 



close systems, but up to half of them could be binaries of 
period shorter than 50 days. We cannot obtain a better esti- 
mate, but this result warns that these so-far unstudied sys- 
tems could play an important role, in particular in clusters 
with a strong lack of close EHB binaries (e.g. NGC 6752). 

- The Gl gap separating the EBTl and the EBT2 in the 
CMD appears as a discontinuity for both the close and the 
intermediate-period binary fraction, because the probabil- 
ity that either f c or fi p is the same in the two populations 
is negligible. Nevertheless, a smooth increase in binarity 
with temperature, rather than a discontinuity, cannot be ex- 
cluded by our observations. If the observed higher binary 
frequency among ch emically unpolluted RGB stars in M4 
dD'Oraziet aLl feolO) should be confirmed as a general be- 
havior of multi-populations in GCs, this could represent a 
problem for the frequently proposed scenario in which the 
EBT2 stars are the progeny of the most He-enriched stellar 
generation observed in the MS. 

- A value of f c as high as that of field stars can be excluded, 
but close EHB systems are surely much more frequent than 
in NGC 6752. In contrast, for NGC 2808 we find results very 
similar to NGC5986. Both these clusters are 2-3 Gyr younger 
than NGC 6752, which indicates that an / c -age is present in 
GCs. 

- The similar / c found in NGC 2808 and NGC 5986 also con- 
firms the problems of reconciling the binary scenario with 
the observations in GCs, because the predictions proposed 
so far cannot account for both the extremely low / c of 
NGC 6752, and the much higher values found in younger 
clusters. We argue, however, that a change in the model pa- 
rameters could enable closer agreement to be reached. 
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Table 3. Data of program stars. 



ID 




RA 


dec 


V 


(U-V) 


T e ff 


RV 


Amax 




hh:mm:ss 








K 


km s -1 


km s -1 


15924 


9: 


1 1:41.89 


-64:42:05.7 


17.975 


-0.932 


14600 


- 


- 


54434 


9: 


12:07.09 


-64:50:31.4 


17.890 


-0.371 


14300 


103.0±0.4 


1.18 


9715 


9: 


11:59.57 


-64:47:44.3 


18.772 


-0.497 


18100 


97.2±1.6 


5.15 


53800 


9: 


11:57.47 


-64:50:41.2 


17.604 


-0.352 


13200 


94.0+0.5 


2.76 


14923 


9: 


11:27.83 


-64:46:22.5 


17.471 


-0.260 


12700 


85.7±1.0 


2.24 


7700 


9: 


11:52.41 


-64:49:18.1 


19.168 


-0.510 


20500 


108.3±2.4 


5.80 


13405 


9: 


11:41.82 


-64:49:16.1 


17.807 


-0.289 


14000 


95.7±1.0 


2.53 


7227 


9: 


12:06.11 


-64:49:32.6 


18.172 


-0.403 


15400 


94.1±0.8 


1.75 


8607 


9: 


12-12 06 


-64:48:42.8 


17.545 


-0.254 


13000 


101.9±0.7 


2.02 


57252 


9: 


12:09.32 


-64:49:43.1 


17.450 


-0.219 


12700 


93.7±0.6 


1.61 


53504 


9: 


12:03.67 


-64:50:45.8 


18.348 


-0.478 


16100 


95.5+1.1 


2.43 


10078 


9: 


12:17.15 


-64:47:17.1 


18.894 


-0.573 


18800 


88.5+2.1 


1.80 


11623 


9: 


12:59.52 


-64:42:51.5 


18.997 


-0.659 


19400 


- 


- 


9121 


9: 


12:21.71 


-64:48:18.3 


18.296 


-0.491 


15900 


91.7+1.2 


2.30 


10620 


9: 


12-29 08 


-64:46:20.8 


17.956 


-0.369 


14500 


92.0+0.9 


2.87 


10719 


9: 


12:24.79 


-64:46:09.9 


17.795 


-0.330 


13900 


101.2+1.1 


1.48 


51625 


9: 


12-08 02 

1 Z. . V./0 .V/Z. 


-64:51:15.8 


17.750 


-0.439 


13800 


100.2+0.5 


2.84 


51455 


9: 


11:58.74 


-64:51:18.5 


17.422 


-0.294 


12600 


90.1+1.4 


1.96 


49247 


9: 


12-08 85 


-64:51:53.0 


17.519 


-0.243 


12900 


- 


- 


50965 


9: 


12:24.44 


-64:51:26.1 


17.714 


-0.288 


13600 


81.3+1.4 


1.96 


47936 


9: 


12:07.89 


-64:52:13.3 


17.902 


-0.465 


14400 


93.3+0.5 


1.57 


47593 


9: 


12:09.69 


-64:52:18.3 


17.632 


-0.312 


13300 


93.0+1.0 


0.91 


53681 


9: 


12-22 70 

1 Z..Z.Z.. / v/ 


-64:50:42.8 


19.095 


-0.620 


20100 


89.8+1.8 


1.82 


53679 


9: 


12:11.84 


-64:50:43.0 


17.470 


-0.178 


12700 


92.1+0.7 


1.71 


55759 


9: 


12:32.23 


-64:50:08.8 


19.064 


-0.628 


20000 


79.2+2.8 


4.04 


52126 


9: 


12:11.18 


-64:51:07.6 


18.098 


-0.386 


15100 


89.4+0.9 


2.00 


55692 


9: 


12:15.38 


-64:50:10.2 


17.804 


-0.297 


13900 


98.4+0.6 


2.55 


9519 


9: 


12:29.66 


-64:47:56.3 


19.271 


-0.619 


21200 


99.2+1.9 


7.95 


47443 


9: 


12:33.43 


-64:52:20.4 


17.751 


-0.374 


13800 


100.1+0.9 


2.46 


48326 


9: 


12:23.85 


-64:52:07.1 


18.846 


-0.673 


18600 


95.4+1.4 


2.68 


50088 


9: 


12:37.34 


-64:51:39.8 


17.795 


-0.348 


13900 


95.9+0.8 


2.02 


50410 


9: 


12:28.99 


-64:51:35.0 


18.232 


-0.453 


15600 


- 


- 


52303 


9: 


12-53 68 


-64:51:03.8 


18.872 


-0.532 


18700 


94.5+1.5 


2.81 


51077 


9: 


12:14.46 


-64:51:24.4 


18.978 


-0.740 


19300 


102.8+1.3 


1.46 


46155 


9: 


12:18.44 


-64:52:42.3 


18.169 


-0.506 


15400 


94.7+1.0 


2.13 


44537 


9: 


12:20.59 


-64:53:09.6 


19.249 


-0.661 


21100 


102.8+2.2 


1.88 


41203 


9: 


12:31.09 


-64:54:14.5 


18.063 


-0.398 


14900 


99.6+1.0 


1.40 


46960 


9: 


12:29.17 


-64:52:28.7 


17.840 


-0.364 


14100 


94.7+1.2 


1.13 


46398 


9: 


12-22 04 


-64:52:38.1 


17.585 


-0.695 


13200 


86.3+1.5 


2.03 


50078 


9: 


12:20.46 


-64:51:40.3 


18.014 


-0.368 


14700 


- 


- 


43422 


9: 


12-06 35 


-64:53:29.6 


18.028 


-0.452 


14800 


98.1+1.9 


0.97 


45166 


9: 


12:17.45 


-64:52:59.0 


17.626 


-0.422 


13300 


82.8+1.5 


2.26 


46460 


9: 


12:11.77 


-64:52:37.3 


18.029 


-0.349 


14800 


98.5+0.7 


1.54 


45677 


9: 


12-02 04 


-64:52:50.4 


17.582 


-0.386 


13200 


99.4+0.7 


1.12 


46512 


9: 


12 08 63 


-64:52:36.6 


17.690 


-0.373 


13500 


92.9+0.8 


0.95 


37288 


9: 


12:34.35 


-64:57:33.5 


17.691 


-0.394 


13500 


92.1+0.8 


2.65 


40263 


9: 


12:07.13 


-64:54:41.3 


19.395 


-0.735 


22100 


- 


- 


37744 


9: 


11:58.76 


-64:56:45.5 


19.229 


-0.542 


21100 


100.2+1.7 


1.88 


37345 


9: 


12:12.46 


-64:57:27.0 


17.916 


0.068 


14400 


-150 


- 


39797 


9: 


12:02.88 


-64:54:55.8 


17.778 


-0.381 


13900 


87.6+0.8 


1.57 


39922 


9: 


12:10.73 


-64:54:51.2 


17.412 


-0.228 


12500 


86.2+0.5 


1.67 


58322 


9: 


12:34.08 


-65:01:05.8 


17.683 


-0.413 


13500 


92.3+0.7 


2.29 


47336 


9: 


12:04.59 


-64:52:22.7 


17.782 


-0.237 


13900 


88.5+0.8 


0.87 


42482 


9: 


12:05.53 


-64:53:47.3 


18.712 


-0.893 


17900 


93.6+1.5 


1.60 


44070 


9: 


12:00.17 


-64:53:18.0 


18.567 


-0.395 


17100 


98.0+1.1 


1.79 


45042 


9: 


11:58.06 


-64:53:01.3 


17.969 


-0.418 


14600 


89.5+0.8 


2.73 


43840 


9: 


12:03.54 


-64:53:21.9 


17.698 


-0.386 


13600 


97.3+0.5 


2.05 


44295 


9: 


11:55.54 


-64:53:13.9 


17.507 


-0.326 


12900 


93.4+0.5 


2.61 


39744 


9: 


11:46.36 


-64:54:57.9 


17.402 


-0.288 


12500 


92.7+0.8 


2.25 


39433 


9: 


12:01.85 


-64:55:09.6 


17.915 


-0.477 


14400 
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Table 4. Data of program stars. 



ID 




RA 


dec 


V 


(U- V) 


T e ff 


RV 






hh:mm:ss 








K 


km s 


km s ' 


43287 


9 


1 1 


53.32 


-64:53:32.1 


18.913 


-0.464 


18900 


86.1+1.4 


0.60 


49443 


9 


] i 


51.16 


-64:51:50.0 


17.460 


-0.266 


12700 


87.1+1.1 


1.60 


47972 


9 


] i 


54.83 


-64:52:12.8 


17.514 


-0.262 


12900 


97.5±0.6 


2.34 


30033 


9 


1 1 


33.68 


-64:56:24.9 


19.269 


-0.794 


21200 






43170 


9 


1 1 


48.45 


-64:53:34.1 


17.649 


-0.333 


13400 


88.0±0.5 


1.92 


47037 


9 


1 1 


54.86 


-64:52:27.7 


17.389 


-0.263 


12400 


90.1±0.6 


2.27 


33067 


9 


] i 


42.45 


-64:52:23.4 


17.740 


-0.325 


13700 


96.5±0.8 


1.98 


32516 


9 


1 1 


16.93 


-64:52:53.6 


17.508 


-0.262 


12900 


92.6±1.0 


2.08 


32670 


9 


] i 


24.84 


-64:52:45.5 


18.509 


-0.730 


16800 


88.2±2.0 


4.74 


45980 


9 


] i 


56.65 


-64:52:45.4 


18.354 


-0.363 


16100 


90.1±0.9 


2.28 


49321 


9 


] i 


54.05 


-64:51:52.0 


18.029 


-0.498 


14800 


85.5±0.6 


2.37 


41115 


9 


1 1 


52.59 


-64:54:17.0 


17.782 


-0.435 


13900 


90.8±0.7 


2.03 


50769 


9 


11 


44.40 


-64:51:29.6 


17.859 


-0.383 


14200 


89.6±0.7 


1.75 


35362 


9 


11 


19.47 


-64:50:04.2 


17.882 


-0.243 


14200 


92.9+1.1 


1.59 


34446 


9 


11 


31.00 


-64:51:04.4 


18.962 


-0.539 


19200 


91.4+2.1 


1.82 


32888 


9 


11 


18.78 


-64:52:32.7 


18.315 


-0.370 


15900 






54539 


9 


11 


50.81 


-64:50:29.9 


18.906 


-0.885 


18900 






54675 


9 


11 


46.80 


-64:50:27.4 


17.516 


-0.345 


12900 


84.3+1.0 


1.86 


55009 


9 


11 


52.90 


-64:50:22.1 


17.807 


-0.325 


14000 


86.0+0.9 


0.83 


13467 


9 


11 


22.19 


-64:49:11.5 


18.996 


-0.558 


19400 






53751 


9 


11 


45.60 


-64:50:41.9 


17.597 


-0.268 


13200 


88.2+1.2 


1.18 


51656 


9 


11 


48.40 


-64:51:15.3 


17.754 


-0.349 


13800 


91.0+0.7 


1.55 


52012 


9 


11 


51.63 


-64:51:09.6 


17.414 


-0.295 


12600 


86.8+0.8 


1.06 



